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Introduction. The development of robotics 1n many advanced countries has raised various industries to a high level. 
The demand for robots increases the share of their use in production tasks, mainly in the motor-vehicle and electronics 
industries. Advanced robotics can increase productivity in many industries by 30%, while reducing labor costs. 
Automation of technological processes of electronics production has a positive impact on the use of robots. Robots are 
used in construction, logistics, oil and gas, aerospace, plant engineering and construction, mining, healthcare, etc. The 
authors consider robots from the point of view of their application in medicine for the rehabilitation of musculoskeletal 
patients. This paper describes the mechanisms for controlling the feet and the center of mass of a humanoid robot. 
Materials and Methods. The authors chose the simplest algorithm for searching for the law of motion control of a 
humanoid robot. The robot movement was presented as a reverse pendulum. Using the large kinematic redundancy of 
walking robots, we have developed a way to control the robot in such a way as to bring the dynamics of its movement 
to the reverse pendulum as close as possible. At the same time, the problem of determining the generalized coordinates 
is considered, at which a given position and orientation of the transferred foot and a given position of the projection of 
the center of mass (CM) of the robot onto the reference surface are provided. 

Results. The authors have developed a digitalized automatic control scheme for the movement of the feet and the center 
of mass of the human exoskeleton, which will largely reduce the load on a sick person. 

Discussion and Conclusions. When discussing the results, comparing the data of the tables obtained during the 
calculation, the following conclusion was made. The scheme for controlling the feet movement of a human exoskeleton 
developed by the authors is most effective when designing an automatic scheme for controlling the movement of the 
feet and the center of mass of a human exoskeleton using digital technology, which will largely reduce the load on a 
sick person. 
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Introduction. In many advanced countries, the successful introduction of robotics has raised various industries 
to a high level. The demand for robots has increased the share of their use in production tasks to 25-45 %, mainly in the 
motor-vehicle and electronics industries. Advanced robotics can increase productivity in many industries by 30 %, 
while reducing labor costs by 18—33 %. Modern robotics can significantly change the entire value chain of products. It 
is estimated that there are about 1.8 million industrial robots operating in the world production systems today, 
representing a global market of about $35 billion: the possibilities of robotics continue to grow, and the costs of 
manufacturing robots continue to fall (they have decreased by about 25 % over the past decade). In production, the 
largest number of robots are used for packaging, grabbing and moving (in Russia, almost 40 % of the 1.7 million), and 
this application has the highest annual growth rate (on average, 11 % per year for 2010-2014). The second common 
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application is in the production of cars, where robots are used primarily for welding. The use of robots for assembly is 
also a fast-growing segment (the average annual growth rate in 2010—2014 was 10 %) due to the growing number of 
electronics/electrical industry products that seek miniaturization and require increased accuracy in manufacturing. 

At Russian enterprises, the density of production robotization is more than 20 times lower than the global 
average. According to the statistics of the International Federation of Robotics, in Kazakhstan, as well as in Russia, 
there are only three industrial robots per 10000 workers, while on average there are 69 worldwide, and more than 100 — 
in the leading countries in the field of digitalization. For example, the share of the Russian market of industrial robots is 
only 0.25 % of the global volume, the main consumers are China (27 %), South Korea (15 %), Japan (14 %) and North 
America (about 14 %). There is also a lag in the share of machines with numerical control: in Japan it is more than 
90 %, in Germany and the USA — more than 70 %, in China — about 30 % [1-3]. 

Modern robotics has been successfully developed in Kazakhstan. Modern factories are being built, where 
robots perform monotonous work. The growth of investments in robotics is increasing every year in Kazakhstan. 
Currently, Kazakhstani enterprises are given a chance to reduce the gap with world leaders. The great flexibility and 
intelligence of robots allow them to be used in various industries where they have not been traditionally used, including 
the production of food and beverages, consumer goods and pharmaceuticals. 

New concepts have emerged, such as wearable electronics and virtual reality, which can be used for quality 
control, work instructions, training, workflow management, various operations, security, logistics, and maintenance. In 
addition, the increased accuracy of these technologies provides improving the profitability of the enterprise by 
increasing productivity. And the instability of quality can be reduced by shortening downtime, defects and waste while 
decreasing the lead time [2]. The development of new technologies will significantly change production processes: it 1s 
most effective in those industries where it is important to adapt the product to customer requirements, and components 
are produced in small volumes and have a high cost. Therefore, the production of consumer goods and the motor- 
vehicle industry, medical and aerospace industries are priorities in the Republic of Kazakhstan, especially in the use of 
3D printing technologies, the production of highly reliable medical devices: hearing aids, dental and other prostheses. 

Materials and Methods. Among all types of robotics, the authors distinguish the medical field. They are 
engaged in the development of exoskeletons of bipedal walking machines (BWM) [4-8]. For BWM, it is generally 
accepted to distinguish two types of walking: static and dynamic. At the static walking, the robot movements are so 
slow that it is possible to neglect the forces of inertia. This allows using a control algorithm built on the basis of only 
equations describing the kinematics of the robot. However, the device speed depends on its size, dynamic parameters, 
kinematic scheme. For many walking machines, this is the main method of control. Figure 1 shows an image of an 
exoskeleton that is being developed by the authors. 





Fig. 1. Functional scheme of the exoskeleton: 
1 — backpack with a microprocessor; 2 — arm; 3 — crutches; 4 — trunk; 5 — foot; 6 — leg; 7 — reference surface 


The equations describing the dynamics of the robot are complex and require the use of well-known methods 
for solving variational problems of large dimension. A simpler algorithm for finding the law of movement control of a 
human robot is needed. The authors from the Bauman Moscow State Technical University proposed the simplest way — 
to imagine the robot as a reverse pendulum and, using the large kinematic redundancy of BWM, control the robot in 
such a way as to bring the dynamics of its movement to the reverse pendulum as close as possible [4—5]. 

Let us consider the problem of determining the generalized coordinates at which a given position and 
orientation of the transferred foot and a given position of the projection of the center of mass (CM) of the robot on the 
reference surface 7 are provided (Fig. 1). The position and orientation of the supporting foot will be considered known. 
Let S be a vector of a given position of the robot foot 5 and CM 6, having dimension 8, 

S = (Xp Ys Zp Op, Bre Ye Xcus Vow - 

It is required to determine the vector gq with dimension 14 BWM, which has 14 controlled mobility stages. 

The task is complicated by the fact that BWM has 14 degrees of mobility, the vector S is 8, and there are kinematic 
restrictions in the mobility stages, the external connections imposed on the feet of BWM change during each step. 
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In this case, various methods are used to solve the inverse problem. For example, an iterative method based on 
minimizing the objective function, and a method for solving the inverse problem in increments based on using the 
Jacobi matrix. However, to control the robot, whose kinematics is represented as a kinematic tree not attached to the 
rack, it is necessary to solve the inverse problem for the supporting foot, the transferred foot, and the projection of the 
center of mass on the reference surface [11-13]. 

The objective function used has the form: f=fp;+fr2tfcutfp where fp; and fr. — the components of the objective 
function that determine the feet positions; fy — a component of the objective function that determines the position of 
the projection of the center of mass; fp — a penalty function that allows bringing the solution closer to the optimal one 
by some criterion. 

On the trajectories of the feet movement and the projection of the center of mass, points are selected so that 
they can then be applied to restore the original trajectory using interpolation with a given accuracy. 7; 1s the matrix of 
the actual position of the foot, and 7,” — the specified position matrix. These matrices are equal if any three points that 
do not match in the connected system have, respectively, the same coordinates in the absolute system [5—9]. 


T,'x, =T?'x, 
T,'x, =T)'x, 
T; Po = T X; 
Setting the vectors ‘x1, 'X>, x3 aS: 
x, =r 
‘X, =H 
17 


i,j,k — the directing vectors of the coordinate axes; 
r iS a parameter that determines the ratio of the accuracy of solving the inverse problem in angular and linear 
coordinates 
AT =1,-T° 
we get 
i. = tr (ATHAT ) 


where the matrix H has the form: 


r 0 0 fr 

0 r Or 
H-= 

0 0 r -r 


rer sr 3 


consider the component of the objective function that determines the position of the projection of the center of mass of 
BWM on the reference surface. The coordinates of the projection of the center of mass are determined by the formula: 


Xo 1 — 0 
S= a > m,PT, T oxai 4 
en M i=6 


where M — robot weight; N — the total number of degrees of mobility of the human robot; m — mass of an individual 
link; roy — vector of the CM link; P — the projecting matrix. 


1 0 0 0 
O 1 0 O 


The index 7 changes from 6, because the links of the kinematic tree with smaller numbers are fictitious and have zero 








mass. If AS, = (X.., Yen y is the target position of the projection of the center of mass, and AS, =S.,-S., , 
then f.,, =AS* 


CM 


Penalty function f, = Ag’ AAg, where g — change of the generalized coordinates when moving to a new point; A — a 


diagonal matrix of weight coefficients. 

Physically, the penalty function fp is proportional to the work performed by the drives when moving to a new 
point [10]. 

When considering the robot movement as the motion of an inverse pendulum, the following equations were 
adopted (Fig. 2): 


_ Xo 
x-x, =lsina, a= hear 


aa a ae LE: 
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Fig. 2. The motion of the reverse pendulum with a constant height of the center of mass 


The pendulum motion along the horizontal x axis is described by the equations 





Poe 
mx = mg 
8 y= 
h h 
take X 9 = 0 and denote 
eae 
l 


We obtain a homogeneous linear differential equation 
x(t)=ce™ +c,e 


t 


Substituting the initial conditions, we get 


a AX. — Xp 
= ane, 
2re 
_ AX, +X, 
C, Xo 
2re 


The equation will take the form 
x(t) = AX, — Xo oe Mt) 4 AX, + Xp eiltto) 
2. 2 

At small angles of inclination of the pendulum, the solution to this equation is close to the solution to the nonlinear 
one written for a pendulum of constant length. Under different conditions, only the expression A [3, 10, 11, 12] will change. 

Research Results. When analyzing the calculations with various options, it was concluded that the matrix 
method is quite time-consuming, takes a lot of time, and does not give measurement accuracy, whereas modern 
methods provide considering the dynamics of the foot movement and the center of mass with greater accuracy, and 
divide the study not into eight or twelve points, but into more, bringing closer to more accurate results. 

Discussion and Conclusions. When discussing the results, comparing the data of the tables obtained during 
the calculation, it was concluded that the scheme for controlling the feet movement of a exoskeleton-human developed 
by the authors is most effective when designing an automatic scheme for controlling the movement of the feet and the 
center of mass of an exoskeleton-human using digital technology. As a result of the obtained equations, an automatic 
control scheme using digital technology was developed, which is shown in Fig. 3 and works as follows. 
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Fig. 3. Block diagram of the exoskeleton control 
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The control device (CD) sets the measurement cycle of the T7- and generates a pulse of the first cycle with 
duration of 77. During the time interval 7;, the electronic key SA is in the initial state, and the measured position of the 
feet and CM through the input device JD goes to the integrator and informs the capacitor (C) the amount of electricity 

1 
GU: 

At the end of the pulse 77, the CD opens the SD key, and the counter Co starts receiving pulses from the 
generator G. At the same time, the SD key is transferred to the second stable state, and the reference value A of the 
reverse polarity is received by the integrator. The capacitor C is discharged to the initial state during the time interval 


Xx 


T 
1360: era 


At the end of the discharge of the capacitor, the voltage at both inputs SS are equal to zero, the SS gives the 
command to open S. The receipt of pulses to the counter stops. Their number determines the 7,. Since the amount of 
electricity during the charge and discharge of the capacitor is the same, then 


T 
U, = re 


1 
Therefore, the measurement result is proportional to the parameters. The time interval 7, does not depend on the time 
constant of the integrator, i.e., no chains with highly stable elements are required to implement the double integration 
method. The duration of the first cycle of integration 7; and the value of the reference values can be kept constant with 
high accuracy, and therefore the error of converting voltage into a time interval is insignificant with this method. The noise 
immunity at the multiplicity of the time of the first integration cycle T; to the interference period reaches 60 dB or more. 
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